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Abstract 

We analyze the peculiar velocity field for 2400 fiat spiral galaxies selected from an infrared sky 
survey (2MFGC). The distances to the galaxies have been determined from the Tully-Fisher 
relation in the photometric J band with a dispersion of 0.”^45. The bulk motion of this sample 
relative to the cosmic microwave background (3K) frame has an amplitude of 199±37 km s“^ 
in the direction 1 = 290° ± 11°, 6 = +1° ± 9°. The amplitude of the dipole motion tends to 
decrease with distance in accordance with the expected convergence of bulk flows in the 3K 
frame. We believe that external massive attractors similar to the Shapley cluster concentration 
are responsible for ~60% of the local flow velocity in the = 0.03 volume. 
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1 INTRODUCTION 

Studying the bulk flows as the deviations from uniform Hubble expansion produced by large- 
scale gravitational effects is of particular interest in testing various cosmological models (Peebles 
1980). The peculiar velocity of a galaxy is a measure of these deviations; apart from the Hub¬ 
ble distance determination, an independent distance determination is required to determine 
this velocity. The Tully-Fisher (TF) relation (Tully and Fisher 1977) and the Faber-Jackson 
(FJ) relation (Faber and Jackson 1976) together with its improved version, the Fundamental 
Plane (FP) method, remain the main popular methods for determining the distances to spiral 
and elliptical galaxies, respectively. The accuracy of determining individual distances by these 
methods is ~20%, and the distances at which they are applicable can reach several hundred 
Mpc. In recent years, more accurate distance determination methods where the absolute mag¬ 
nitudes of type-la supernovae (SN la) (Riess et ah 1997) or the surface-brightness fluctuations 
(SBF) of galaxies (Tonry et ah 2000) are the distance indicators have been developed. Dis¬ 
tances of several hundred Mpc can be estimated from SN la to within 8% (unfortunately, the 
number of measured supernovae is small). Distances as large as several tens of Mpc can be 
estimated by the SBF method with an accuracy of 10-20%. 

A huge body of observational data (see, e.g., Willick (2000) Courteau and Dekel (2001), 
Zaroubi (2002) and references therein) pertaining separately to E and S galaxies, field galaxies. 
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cluster members, etc. has been accumulated over 20 years of active studies of non-Hubble 
flows. The works of different groups also differ by the depth of the samples studied and the 
peculiarities of the sky distribution of galaxies. The magnitude of the velocity and the direction 
of the bulk motion can be determined by measuring the radial peculiar velocities of galaxies. 
For the standard ACDM model, the bulk velocity in the cosmic microwave background (3K) 
frame is expected to approach zero as the volume under consideration increases. Therefore, 
measuring the dipole component on various scales is needed to determine the volume in which 
the flow converges. At present, there is no contradiction between the velocities (~ 220 km s“^ 
) and the directions of the bulk motions of galaxies estimated by different anthors on a scale 
of ~60 Mpc. At this depth, the data for E and S galaxies, held galaxies, and clnsters are in 
good agreement, within the error limits, irrespective of the distance determination method (see 
Table 1 in the review by Zaronbi (2002)). 

Two gronps of results have been obtained at distances of 100-150 Mpc. Using 522 spiral 
galaxies in Abell clusters. Dale et al. (1999) determined the magnitnde of the bnlk velocity, 
V = 75±92 km s“^, and the apex position, I = 282° and b = 25°(r ~ 150 Mpc), by the 
TF method. Using 85 SN la, Riess et al. (1997) derived the apex position, I = 289° and 
b = —2° , at a magnitnde of the velocity that differs only slightly from zero. The results 
obtained by Colless et al. (2001) using the FP method apply to distant E galaxies in clusters 
located at distances of Hor = 6000 — 15000 km s“^ (EFAR) toward Hercules-Corona Borealis 
and Persens-Pisces-Cetns. The anthors fonnd the bnlk velocity in the volnmes stndied to be 
consistent with the hypothesis abont its zero valne at a 5% significance level. The above results 
are in good agreement: on scales of 100-150 Mpc, the bnlk velocity is 0-200 km s“^ and the 
apex positions coincide, within the error limits, and are close to the direction of the excess of 
point sources in the sky from the IRAS PSCz catalog (Sannders et al. 2000). It is important 
that the convergence of bnlk flows in this gronp of results was obtained for different types of 
objects and using different distance determination methods. 

The second gronp of results is characterized by a high magnitude of the bulk velocity and a 
signihcant spread in apex positions. Based on photometry for the brightest galaxies in clusters, 
Lauer and Postman (1994) determined the parameters of the motion relative to the 3K frame 
for 119 Abell clusters of galaxies with velocities up to 15 000 km s“^ (LP): V = 689±178 km s“^ 

, I = 343°, b = 52° ± 23° . Willick (1999) found the parameters of the bnlk motion for 15 Abell 
clusters of galaxies located at a distance of 120 Mpc, V = 720 ± 280 km s“\ I = 272°, b = 10° 
(using the TF method). Hndson et al. (2004) fonnd the parameters of the bnlk motion in the 
3K frame for 56 Abell clusters (the SMAC sample) up to Hqt = 6300 km s“^ with an effective 
sample depth of = 6300 km s-1 : V = 687 ± 203 km s“^, / = 260° ± 13°, 5 = 0° ± 11° . 
The anthors believe that, apart from the Shapley Concentration, even more distant attractors 
must provide such a high velocity. Hudson (2003) attempted to reconcile the conflicting data 
at depths from 6600 to 11000 km s“^ by assnming that the sample sparseness and small size 
lead to large errors in the peculiar velocities. If the LP sample is excluded, then the results 
with high and low bnlk velocities can be reconciled at a 2a level. The sample that combines the 
SMAC survey with other surveys yields a bulk flow with a velocity of 350±80 km s“^ toward 
I = 288° ,6 = 8° (Hudson 2003). 

What the converge scale is remains an open qnestion nntil new observational data based on 
mnch more complete and homogeneous samples are obtained. There are advantages and dis¬ 
advantages in nsing clusters of galaxies to study the bulk motions. The distance determination 
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for a cluster is based on averaging the distances to its individual members, which reduces the 
error. However, the number of measured galaxies in clusters is generally small. In addition, 
including the cluster held galaxies remains a possibility. The number of clusters measured by 
various authors ranges from 20 to 100. Even if the cluster samples are distributed more or less 
uniformly over the sky, they are very sparse due to the small size. The measurement errors 
of the peculiar velocities for held galaxies are fairly large and increase with distance and sam¬ 
ple incompleteness. This shortcoming is compensated for by the possibility of increasing the 
sample to hundreds or even thousands of galaxies. The Mark III catalog (Dekel et ah 1999) is 
an example of combining several catalogs of peculiar velocities to determine the parameters of 
the bulk motion. This catalog includes ~3000 E and S galaxies located at i? < 60 Mpc. The 
Mark III galaxies move relative to the 3K frame with a bulk velocity oiV = 370 km s“^ toward 
/ = 305°, b = 14°. This velocity seems fairly high compared to V obtained for other samples 
located in the same volume. The authors believe that an uncertainty in the mutual calibration 
of the heterogeneous samples may be responsible for the discrepancy. 

Thus, an extensive deep homogeneous sample of galaxies uniformly and densely distributed 
over the sky is needed for a more dehnite solution of the question regarding the convergence 
scale of bulk flows. The edge-on flat galaxy catalog, RFGC (Karachentsev et al. 1999), satishes 
these conditions. The RFGC includes 4236 galaxies with optical angular diameters a > 0.6' 
and apparent axial ratios a/h > 7. The RFGC is currently the most uniform all-sky catalog of 
edge-on late-type spiral galaxies. About 3000 objects in the catalog have fairly reliable J, H, K 
photometry in the 2MASS survey (Cutri and Skrutski 1998), and ~1200 objects have estimates 
of the radial velocities 14 and the hydrogen line widths We constructed the 2MASS-TF 
relations for a sample of 1141 galaxies and, after cleaning, determined the parameters of the 
bulk motion for galaxies with V^k < 12000 km s“^ using 921 galaxies: V = 226 ± 62 km s“^, 

1 = 295° ± 16°, and b =-2° ± 13° (Kudrya et al. 2003). 

Supplementing the 2MASS-RFGC sample to study the bulk flows suggests performing new 
measurements of 14 and H4o for RFGC galaxies. In particular, this program has been imple¬ 
mented with the 100-m Effelsberg radio telescope by Huchtmeier et al. (2005) and Mitronova 
et al. (2005) since 2001. Compiling a deeper homogeneous sample of late spirals based on the 
2MASS catalog (2MFGC) is a fundamentally different possibility (Mitronova et al. 2004). 

2 THE 2MFGC CATALOG 

The successful application of the infrared (IR) Tully-Fisher relation to RFGC flat galaxies to 
determine the parameters of the bulk motion (Karachentsev et al. 2002; Kudrya et al. 2003) 
was a justihcation for compiling the new 2MFGC catalog. The 2MASS survey is known to have 
a low sensitivity with respect to late-type galaxies (Jarrett 2000), since the periphery of the 
disks of spiral galaxies is unseen on isophotes fainter than Kg = 20™. 

Studying the near-IR properties of RFGC galaxies and comparing them with optical data 
allowed the selection criteria to be worked out when compiling the new flat galaxy catalog based 
on 2MASS. Thus, for example, comparison of the IR and optical diameters of RFGC galaxies 
showed that the 2MASS diameters are, on average, half the standard optical radii. For RFGC 
galaxies, the optical axial ratios a/b cover the range (7-21) with a median of 8.6, while the 
corresponding IR axial ratios he within the range (1-10) with a median of 4.1 (Karachentsev 
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et al. 2002). 

The principles of selecting galaxies for the new catalog from the original 2MASS-XSC 
catalog were described in detail by Mitronova et ah (2004). The main parameter in the selection 
was the apparent axial ratio given in the XSC (Extended Sources Catalog). The all-sky 2MFGC 
disk galaxy catalog (Mitronova et al. 2004) contains IR photometry and identihcations in the 
LEDA and NED databases for 18020 objects with axial ratios a/h > 3 in the XSC-2MASS. 
Analysis of the distribution of galaxies in axial ratio in the new catalog yields a median value 
of a/h = 4, which is close to the median value of a/h = 4.1 for the RFGC galaxies identihed 
in the XSC-2MASS. The two catalogs are similar in morphological composition and contain 
~80% of the late-type spiral galaxies. 

Below, we use the following photometric parameters from the 2MFGC: J/e , the elliptical 
Kron magnitude; r^e , the elliptical Kron radius (in arcseconds) in the K band; a/6, the axial 
ratio averaged over all J , FT , and K bands; and IC, the concentration index (the ratio of the 
radii within which 3/4 and 1/4 of the light from the galaxy is concentrated). 

3 A SAMPLE FOR DETERMINING THE PARAME¬ 
TERS OF THE BULK MOTION 

Compiling a sample of 2MFGC galaxies with estimates of the heliocentric radial velocity V/ 
, the HI line width W^q , W 20 , or the maximum rotational velocity Vm are the subject of a 
separate paper. The radial velocities are known for 5536 of the 18 020 2MFGC galaxies, and 
estimates of HAo , W 20 , or Vm are also available for 2765 objects. The widths HAo (N=2276) and 
HAo (N=1825) are represented most completely in the sample. Both HAo and W 20 estimates are 
available for 1740 sample galaxies. Optical measurements of Vm are available for 445 galaxies; 
32 of them also have W 50 estimates. To choose the basic width, HAo or W 20 , in determining the 
parameters of the bulk motion, we hrst constructed the standard twoparameter TF relations 
for both widths. We take the TF relation in the form 

M = Gi+ C'2logIW. (1) 

We calculate the absolute magnitude from the apparent magnitude J/e in the standard way, 

M = J/e-25-5 log r, (2) 

by hrst estimating the photometric distance r [Mpc] using the post-Hubble relation 

r = Vsk/Ho{ 1 — {qo — 1)V3 k/‘^c}, (3) 

where the radial velocity V^k in the CMB frame is calculated from the heliocentric radial 
velocities V/ = cz using the parameters of the solar motion relative to the background from 
Kogut et al. (1993). We take the Hubble constant Hq = 75 km s“^ Mpc“^ and the deceleration 
parameter go = —0.55 in accordance with the cold dark matter (GDM) cosmological model and 
the cosmological constant (Dm = 0.3, Da = 0.7). 

We correct the widths for the cosmological expansion, 

W'^ = W/{l + z), (4) 
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by abandoning the correction for tnrbnlence, which was applied in onr previons paper (Kndrya 
et ah 2003), since it tnrned ont to have no effect on the bnlk velocity parameters. We determine 
the coefficients of relation ( 1 ) by least sqnares with eqnal weights. 

We applied mainly the same data cleaning techniqne that was nsed previously (Kndrya et 
al. 2003) by excluding galaxies with deviations of more than 3 (Ttf in the TF diagram and 
with individual peculiar velocities in the 3K frame V^ec > 3000 km s“^ . In contrast to our 
previous paper (Kndrya et al. (2003), we did not exclude distance galaxies {V^k > 18000 km 
s“^, Hqt > 18000 km s“^) since there were only a few such objects. After the cleaning, we 
obtained samples of N = 1604 galaxies with IF 20 and N = 1972 galaxies with IF 50 . The 


Table 1: Parameters of the Tully-Fisher relation for various samples 


N 

Sample 

axF 

C'l 

C '2 

o-y 

1 

IF 20 , N1825 

1.01 

-5.42 ±0.35 

-6.46 ±0.14 

1925 

2 

W 50 , N2276 

0.97 

-7.14 ±0.28 

-5.86 ±0.11 

1927 

3 

IF 20 , A^1604 

0.55 

-1.14 ±0.23 

-8.15 ±0.09 

992 

4 

IF 50 , Ari972 

0.49 

- 2.10 ±0.18 

-7.86 ±0.07 

974 

5 

IT 50 ± Vm, N2Q8A 

0.92 

-6.80 ±0.25 

- 6.01 ± 0.10 

1921 

6 

IT 50 ± Vm, N2333 

0.49 

-2.55 ±0.17 

-7.69 ±0.07 

1015 

7 

W 50 ± Kn ± hF20) N27Q5 

0.94 

-7.11 ±0.25 

-5.88 ±0.10 

1930 

8 

W 50 ± Vm ± IF 20 ) A^ = 2395 

0.49 

-2.68 ±0.17 

-7.64 ±0.07 

1016 


results of our calculations are presented in the hrst two rows of Table 1 for the complete samples 
of N = 1825 galaxies with IF 20 and N = 2276 galaxies with IF 50 and in the third and fourth 
rows for the cleaned samples. 

In Table 1, utf denotes the standard deviation (in magnitudes) in relation (1), Ci and C 2 
are the coefficients in ( 1 ), and ay is the standard deviation (in km s“^ ) in peculiar velocity 
space. In both cases, the scatter in the TF diagram is seen to be smaller for the widths IF 50 . 
Therefore, we take IF 50 as the basic width and then successively add the data with VA and IF 20 
to the sample with known IF 50 . 

To pass from Vm to IF 50 , we used the relation 

IT5o = 2Id„ +(27.4 ±7.8), (5) 

derived from 32 galaxies with measured IF 50 and Vm- In this way, we enlarged the sample by 
408 objects, thereby bringing it to N = 2684. 

The calculated parameters of the TF relation for the complete N = 2684 sample and for the 
N = 2333 sample cleaned using the adopted technique are given in the hfth and sixth rows of 
Table 1. We see that adding the data with Vm slightly reduced the scatter in the TF relation 
for the uncleaned data and left the scatter for the cleaned data unchanged. 

Our sample contains 1740 objects with measured IF 50 and IF 20 . The relation between these 
widths is shown in Fig. 1. 

After excluding the 68 galaxies that deviate from the linear relation by more than 3a, we 
obtain the relation 

IT 50 = IF 20 -(19.8 ±0.5). ( 6 ) 
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Figure 1: Relation between W^o and The squares mark 68 galaxies that deviate from the 
linear relation by more than 3a. 
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Figure 2: Tully-Fisher relations for the initial (a) and cleaned (b) samples. 


The hypothesis that the coefficient of the linear correlation between the widths is equal to 
unity is conhrmed by the Fisher test for an orthogonal regression. Adding the data with 

fF 2 o enlarged the sample by 81 galaxies. Parameters of the TF relation (1) for the complete, 
N = 2765, and cleaned, N = 2395, samples are given in the seventh and eighth rows of Table 1. 

The Tully-Fisher relations for the initial and cleaned samples are shown in Fig. 2. The 
dashed lines in this figure correspond to the linear regression constrncted by least sqnares. 

We give Fig. 2a to illustrate the quality of the initial observational data. The region where 
most of the galaxies are concentrated is clearly seen, bnt an appreciable fraction of them (13%) 
exhibit a large scatter in the diagram. As we see from Fig. 1, this scatter is partly (by abont 
1/3) attribntable to the low qnality of the HI spectra. Comparison of Figs. 2a and 2b shows 
that our data cleaning criteria are efficient. 

In the next section, we determine the parameters of the bnlk motion for the cleaned N = 
2395 sample. 
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Figure 3 presents the sky distribution of 2395 galaxies in Galactic coordinates. The galaxies 
with measured IT 50 (N = 1977), (N = 353), and W 20 (N = 65) are marked by hlled circles, 
squares, and crosses, respectively. As we see, combining all of the data improves appreciably 
the uniformity of the sky distribution. 

4 PARAMETERS OF THE BULK MOTION FOR 2MFGC 
GALAXIES 

4.1 Basic Formulas and General Results 

We used a simple two-parameter TF relation to compile our main sample of galaxies. To 
calculate the parameters of the bulk motion, we generalize relation ( 1 ) by including additional 
regressors. By an exhaustive search for various parameters of galaxies, we found that three 
regressors are most signihcant according to the Fisher test: the radius r/g , the axial ratio a/b, 
and the concentration index IC. This linear relation is taken in the form 

M = Cl -\- C 2 ■ log -\- C *3 ■ logrje + C *4 ■ a/h -\- C^ ■ IC. (7) 

In this case, formulas (2)-(6) remain valid. We calculate the individual peculiar velocity in 
the post Hubble approximation, 

Vpec = VsK - Hor{l + (go - l)i7or/2c}. (8) 

We use the set of peculiar velocities to calculate the orthogonal components of the dipole 
bulk velocity V, 

Vpec,i = V ■ Ci AVi (9) 

by minimizing the sum of the squares of the “noise” peculiar velocity component Vi {i is the 
galaxy number in the sample). Here, e) = (cos/j cos&j, sin/j cos 6 j, sin 6 j) is a unit vector of 
direction of galaxy i in a reference frame associated with the Galactic coordinates I and b. 

We calculated the errors in the magnitude of the bulk velocity V and the apex position /, b 
as follows. First, we determined the diagonal elements Byy , Bu , and of the covariance 
matrix B in the basis {eV, e), e),} and then AH = A/ = arcfg{(Hjj)^/^/H}, Ab = 

arctg{{Bi,bf/‘^/V}. 

The calculated bulk velocity parameters for the N = 2395 sample and for the samples of 
2MFGG galaxies with the additional restrictions V^k < 8000 km s“^ and < 10000 km s“^ 
are given in the hrst three rows of Table 2. Here, I 4 , I 4 denote the orthogonal components 
of the dipole bulk velocity, V is its magnitude (all in km s“^ ), I and b are the Galactic 
coordinates of the apex, and F is the signihcance coefficient according to the Fisher test for 
the bulk velocity vector. Gomparison of the hrst row in Table 2 and the eighth row in Table 1 
shows that including the additional regressors reduces the scatter in the TF diagram by 6%. 

To test the robustness of the results obtained and to illustrate the signihcance of the post- 
Hubble corrections in formulas (3) and (8), we calculated the parameters of the hve-parameter 
regression and the bulk velocity by taking a linear Hubble law. The results for the N = 2395 
sample are given in the fourth row of the table. We see that the apex position changed only 
slightly, while the magnitude of the bulk velocity increased by ~ 5%. Our calculations for other 


7 



90 



- 90 ® 


Figure 3: Distribution of 2395 2MFGC galaxies in Galactic coordinates. Galaxies with mea¬ 
sured VFso , Vm , and hF 2 o are denoted by filled circles, squares, and crosses, respectively. 

samples also confirm that allowance for the relativistic redshif-distance relation slightly reduces 
the estimated magnitude of the bulk velocity. 

We also checked how the sample cleaning based on the five-parameter regression (7) affected 
the parameters of the bulk motion. The galaxy selection criteria were the same as those for 
the two-parameter regression. The results of our calculations for the new N = 2377 sample are 
presented in the hfth row of Table 2. Gomparison of the data in rows 5 and 1 shows a tendency 
for the magnitude of the velocity to decrease for a more rigorous sample cleaning. At the same 


Table 2: Parameters of the dipole bulk velocity 


Sample 

gtf 

cry 

U 



V 

1 

b 

F 

A11,N2395 

0.46 

962 

109±37 

-196±37 

-20±31 

225±36 

299 ± 10 

-5 ±8 

13 

VsK < 10000, N2205 

0.48 

914 

125±37 

-201±37 

-25±30 

238±36 

302 ±9 

-6 ±7 

15 

VsK < 8000,N1943 

0.49 

856 

87±37 

-206±37 

-22±30 

225±36 

295 ± 10 

-6 ±8 

14 

LHL,N2395 

0.46 

970 

110±37 

-206±38 

-28±31 

235±37 

298 ±9 

-7 ±7 

14 

N2377 

0.45 

968 

67±37 

-187±38 

2±31 

199±37 

290 ± 11 

1 ±9 

10 


time, the data in Table 2 indicate that the sought-for values of V, /, and b agree, within the 
error limits. Figure 4 presents the depth dependence of the magnitude of the bulk velocity 
(dots with bars, left vertical scale) and the sample size to a given depth (open diamonds, right 
vertical scale). 

We see that the magnitude of the velocity decreases for distances from 3000 to 8000 km 
s“^, as expected in the standard model for the formation of structures in the Universe with 
the cosmological constant and cold dark matter. The sample size increases roughly linearly to 
V^K = 8000-10000 km s“^ and then is rapidly saturated. After 8000 km s“^, the magnitude of 
the velocity (230-240 km s“^) and the apex position (which is not given here) do not change, 
within the error limits. Thus, our results yield the magnitude and direction of the bulk velocity 
for 2MFGG galaxies to an effective depth of ~ 10 000 km s“^. 
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Figure 4: Magnitude of the bulk velocity (dots with bars, left vertical scale) versus depth and 
the sample size to a given depth (open diamonds, right vertical scale). 


4.2 Parameters of the Bulk Motion for Spherical Shells 

We also calculated the parameters of the bulk motion separately for four subsamples of the 
complete N = 2395 sample by breaking it down into shells: OkVsk < 3000 km s-^ Si; 3000 < 
y^K < 6000 km s“^ , S 2 ; 6000 < V^k < 9000 km s“^, S^; and 9000 < V^k < 20000 km s“\ 5*4 . 
The volumes of the first, second, third, and fourth subsamples include, respectively, the Local 
Supercluster (/ = 284°, b = 74°), the Great Attractor (Hydra, / = 270°, & = 26°; Cen 30+45, 
I = 302° , b = 22°), the Coma supercluster (/ = 58°, 6 = 88°), and the Shapley Concentration 
(/ = 312°, 6 = 31°). 

The sky distribntions of galaxies in the shells in Galactic coordinates are shown in Fig. 
5. The bonndary of the shell (in km s“^ ) and the nnmber of galaxies in it are given in the 
npper left corner of each panel. The regions of enhanced density corresponding to well-known 
snperclusters are clearly distinguished in the individnal panels. Table 3 presents the calculated 
dipole bulk velocity. The parameters of the Tully-Fisher relation were determined for each shell 
separately. The content of the columns is the same as that in Table 3. As follows from these 
mntnally independent data, the apex position and the bnlk velocity amplitude in each shell 
remain approximately the same within {l-2)a. The amplitnde is at a minimum, (157+45) km 
s“^, in shell S 2 , where the Great Attractor is located. We see no evidence of the flow toward 
the Virgo cluster within the snbsample Si with V^k < 3000 km s“^; i.e., the entire volnme of 
the Local Snpercluster moves approximately in the same direction as the remaining objects of 
our sample. 


4.3 Parameters of the Bulk Motion in the Frame of the Local Group 

As we noted above, ~13% of the galaxies in the initial sample deviate signihcantly from the 
linear TF regression (see Fig. 2a). The low qnality of the 2MASS photometry or the HI 
spectra for these galaxies is usnally responsible for the deviations. Occasionally, the deviation 
from the TF relation is caused by the peculiar structure of a galaxy or the tidal perturbation 
from its close neighbors. We hoped to eliminate these cases by applying the three sigma rule 
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and restricting the admissible absolute value of the peculiar velocity by the obvious physical 
limit Vpec < 3000 km s“^. As we see from relations (l)-(3) and (8), the adopted conditions for 
eliminating “bad” galaxies depend on the reference frame in which the distance and absolute 
magnitude of the galaxy are determined from its radial velocity. Therefore, in principle, the 
dipole parameters can contain a systematic error caused by the a priori choice of the 3K frame 
for estimating the distances to galaxies. 

According to Kogut et al. (1993), the Sun moves relative to the 3K frame with a velocity 
of V^K = 369.5 ± 3.0 km s“^ in the direction I = 264.4° ± 0.3°, b = 48.4° ± 0.5°, while according 


Table 3: Parameters of the dipole bulk velocity for galaxies in the spherical shells 


N 

Sample 

(Jtf 

cry 

K 


K 

V 

1 

b 

F 

1 

S'i,N509 

0.57 

438 

236 ± 42 

-276±36 

27±27 

364±38 

311 ±6 

4±4 

31 

2 

S'2,N956 

0.37 

738 

38 ±43 

-145±44 

-47±41 

157±45 

285 ± 16 

-17± 14 

4.4 

3 

S'3,N645 

0.30 

1003 

118 ±73 

-315±83 

-147±61 

367±85 

291 ± 11 

-24 ±9 

6.4 

4 

S'4,N285 

0.25 

1187 

157± 136 

-287±139 

-75±109 

335±143 

299 ± 22 

-13± 18 

1.8 


to Karachentsev and Makarov (1996), the Local Group centroid has a velocity of V^k (LG) 
= 634 ± 12 km s“^ relative to the 3K frame in the direction I = 269° ± 3°, 6 = 28° ± 1° . 
The orthogonal components of this velocity 14,1^, 14 in the Galactic coordinate system are 
(—8, —559, 298) km s“^ . To check the amplitude and direction of the bulk motion for the 
possible presence of a systematic error, we calculated the parameters of the dipole solution for 
2MFGG galaxies from their velocities relative to the Local Group centroid and compared them 
with those obtained above in the 3K frame. 

For N = 2404 galaxies that satished the three sigma rule and the condition 14ec < 3000 
km s“^ in the LG frame, we obtained a standard deviation relative to the two-parameter TF 
relation of axp = 0.49™ or ay = 986 km s“^ . The bulk motion of these galaxies in the LG 
frame has a velocity oiV = 372 ± 33 km s“^ in the direction / = 76° ± 6° , 6 = —40° ± 6° ; the 
orthogonal components of this velocity are . Since 

%k{2MFGC) = Vlg{‘^MFGC) + V^k{LG), 

returning to the 3K frame, we obtain the orthogonal bulk velocity components (-1-60, —284, 
-|-57) or the following new amplitude and direction of the bulk motion: V = 296 km s“^, 
I = 282°, b = 11°. Gomparison with the previous data (see row 1 in Table 2) shows that the 
systematic difference in the parameters of the bulk motion due to different galaxy elimination 
procedures lies within 2a of the random errors. However, it should be noted that the systematic 
difference in Vx,Vy,Vz is more signihcant in the close volumes Si and S2 ■ 

4.4 The Peculiar Velocity Field 

The distribution of the peculiar velocities for 2MFGG galaxies dehned by relations (7) and (8) 
is presented in Galactic coordinates in Fig. 6. Galaxies with positive and negative peculiar 
velocities are indicated by the open and hlled circles, respectively. This discrete distribution was 
averaged by a Gaussian Liter with a 20° window. Equal mean peculiar velocities are indicated 
by the lines at 50 km s“^ steps. The heavy line corresponds to a zero peculiar velocity; negative 
Vpec are indicated by the dashes. 
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As follows from these data, the peculiar velocity distribution of 2MFGC galaxies is asymmet¬ 
ric: the maximum positive isovelocity is -1-350 km s“^, while the maximum negative isovelocity 
is only 150 km s“^. The Vpec distribution shows a distinct dipole pattern. The region of positive 
mean Vpec appears doubly connected, with the primary peak being approximately in the zone 
where the Hydra-Centaurus, Norma, A3627 clusters and the Shapley Concentration of clusters 
are located. The secondary positive peak with an amplitude of 150 km s“^ is identified with the 
Abell cluster A400 (170° , —45° ) and A539 (196°, —18°). The region of negative mean peculiar 
velocities forms three broad shallow troughs with the coordinates of their centers (120°, -|-40°), 
(80°, —30°), (200°, -|-30°). The hrst two of them he not far from the well-known Void in Bootes 
and the Local Void. In general, the relief of the crests and troughs in Fig. 6 agrees well with 
that obtained from a less representative sample (N = 971) of RFGC galaxies with 2MASS pho¬ 
tometry (Kudrya et ah 2003). The reproducibility of general features in the Vpec distribution 
from two essentially independent observational samples gives hope that the distribution of the 
gravitational potential in the local Universe can be reconstructed from the available peculiar 
velocity map. 

5 DISCUSSION AND CONCLUSIONS 

To estimate the representativeness and quality of the galaxy sample used to determine the 
dipole parameters, we introduced the concept of goodness, G = (iV/100)^/^/(TrF , where N is 
the number of galaxies in the sample, and TF is their dispersion (in magnitudes) in the TF 
diagram (Kudrya et ah 2003). Of course, apart from N and the quality of the sample is 
also characterized by the uniformity with which it covers the entire sky. However, allowance 
for the nonuniform distribution of galaxies over the sky requires more complex approaches. At 
present, there are only hve samples whose goodness exceeds G = 5. They are all listed in Table 
4 in order of increasing G. Golumns 1, 2, and 3-5 give the galaxy type, its goodness G, and 
the dipole amplitude and direction in Galactic coordinates, respectively; the last column gives 
a reference to the data source. 

As follows from these data, the sample of 2MFGG galaxies being discussed has the best 
goodness. The parameters of the 2MFGG dipole are in good agreement with their weighted 
mean values for all hve samples: V = 225 ± 47 km s“^, / = 295° ± 3°, & = -|-6° ± 5° . These 
amplitude and direction of the bulk motion can be currently taken as the standard ones for the 
bulk how of galaxies relative to the GMB on a scale of ~100 Mpc. It should be noted that this 
direction of the dipole is almost opposite to the direction of the rotational velocity of the Sun 
relative to the Galactic center, V = 220 km s“^ , / = 90°, 6 = 0°. For this coincidental reason, 
the large-scale how of galaxies had gone unnoticed until the 1980s. 

Recall that we compiled and used catalogs of edge-on spiral galaxies to analyze the bulk 
motions of galaxies. The selection of such galaxies by their optical parameters (FGG) led 
to the following apex parameters: V = 300 ± 75 km s“\ I = 328° ± 20° , and b = 78° ± 158° 
(Karachentsev et ah 1993). An increase in the sample size in combination with the passage from 
optical to IR photometry yielded new apex parameters: V = 199 ± 61 km s“^, I = 301° ± 18°, 
and b = —2° ± 15° (Kudrya et ah 2003). Using the new 2MFGG catalog compiled from the 
2MASS survey provided a more uniform distribution of galaxies over the sky. The amplitude 
and direction of the bulk how calculated here, V = 199 ± 37 km s“^ , I = 290° ±11°, and 
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Figure 6: Peculiar velocity distribution for 2395 2MFGC galaxies in Galactic coordinates. 
Galaxies with positive and negative peculiar velocities are indicated by the open and hlled 
circles, respectively. The velocities of the galaxies were averaged by a Gaussian hlter with a 20 
window. The lines of constant mean peculiar velocity are indicated at 50 km s“^ steps. The 
heavy line corresponds to a zero peculiar velocity; negative Vpec are indicated by the dashes. 

b = +1° ± 9°, are in satisfactory agreement with our previous dipole solutions (Karachentsev 
et ah 2000). Thus, our parameters of the local large-scale flow turned out to be stable with 
respect to the selection of galaxies by their optical or IR parameters, the use of the optical or 
IR luminosities in the TF diagram. The dipole parameters also change little if other regressors 
are included in the standard TF relation between absolute magnitude and 21-cm line width; 
the axial ratio, the surface brightness, the concentration index, and quadratic combinations of 
these regressors. The manner of eliminating galaxies by their peculiar velocity (in the 3K or 
LG frame) introduces a systematic errors within 2a of the random errors. 


Table 4: Goodness for various samples 


Sample 

G 

V 

1 

b 

Reference 

Spirals in clusters 

6.0 

75 

289 

25 

Dale et ah (1999) 

SNIa 

6.8 

206 

290 

0 

Radburn-Smith et al. (2004) 

RFGG 

7.4 

199 

301 

-2 

Kudrya et al. (2003) 

Mark III 

8.6 

370 

305 

14 

Dekel et al. (1999) 

2MFGG 

10.9 

199 

290 

1 

This paper 

Mean 

— 

225 

295 

6 

— 


The apex of the bulk flow is located roughly in the same region as the centroids of IRAS 
sources (/ = 258°, 6 = -|-30°), the centroid of 2MASS galaxies (/ = 278°, 6 = -|-38° ), and the 
concentrations of rich Shapley clusters (/ = 315°, & = -|-30° ). If judged by the pattern of 
decrease in the flow amplitude in Fig. 4 with linear scale of the Local Volume, then objects 
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(attractors) located farther than 50 Mpc could be responsible for ~60% of the amplitude of the 
observed bulk motion of 2MFGC galaxies. 
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